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Abstract 

We study the locations of emission regions through modelling the quasi-simultaneous multi-frequency 
spectral energy distributions of 21 Fermi flat spectrum radio quasars (FSRQs) in the frame of a multi- 
component one-zone leptonic model. In our calculations, we take the detailed broad line region (BLR) 
structure into account and discuss the effect of the uncertainty of the BLR structure on constraining 
the location of the emission regions for each FSQR, meanwhile we also include both the internal and 
external absorptions. Our results indicate that (1) the contribution of external Compton-BLR component 
is important to 7-ray emission, and the energy density of external target photon fields depends on the 
location of the emission region, which can be derived through reproducing the observed 7-ray emission, 
and (2) the emission regions of FSRQs with relative low accretion disk luminosity lie in the region of 
7.9 X lO^s - 1.3 X 10^* cm (300 - 4300 Schwarzschild radh) from central black hole, and for FSRQs with 
high accretion disk luminosity, the emission regions locate in the larger region of 2.6 x 10^'' - 4.2 x 10^* cm 
(300 - 5600 Schwarzschild radii). 

Key words: galaxies: active - radiation mechanisms: nonthermal - quasars: general - 7-rays: obser- 
vations 



1. Introduction 

Blazars are the most extreme class of active galac- 
tic nuclei (AGNs). Their spectral energy distributions 
(SEDs) arc characterized by two distinct bumps: the first 
bump located at low-energy band is dominated by the 
synchrotron emission of relativistic electrons which origi- 
nates in a relativistic jet, and the second bump located at 
high-energy band could be produced by inverse Compton 
(IC) scattering (e.g., Bottcher 2007). Various soft pho- 
ton sources seed the synchrotron self-Compton (SSC) pro- 
cess (e.g., Rees 1967; Maraschi et al. 1992) and exter- 
nal Compton (EC) process (e.g., Dermer & Schlickeiscr 
1993; Sikora ct al. 1994) in the jet to produce 7-rays. 
It should be noted that hadronic models have also been 
proposed to explain the multi-band emissions of blazars 
(e.g., Mannheim 1993; Miicke et al. 2003). At present, 
the detectors including the ground-based and satellite- 
based ones can provide us the high quality data to con- 
struct SEDs. Therefore, the simultaneous and quasi- 
simultaneous SEDs become the basic tools to study the 
physics properties of blazars, such as magnetic field, the 
size of emission region, Doppler factor and so on. 

The Large Area Telescope (LAT) onboard the Fermi 
satellite provides unprecedented sensitivity in the 7- 
ray band. In its two years operations, 987 blazars 
were detected, including 360 flat spectrum radio quasars 
(FSRQs), 423 BL Lacs and 204 unknown type blazars 
(Ackerman et al. 2011). With Fermi, Swift and other ob- 
servatories. Abdo et al. (2010a) assembled high quality 
data of 48 blazars in the first three months LAT sam- 



ple (LAT Bright AGN Sample: LBAS) to build quasi- 
simultaneous SEDs. The data from Swift were col- 
lected in one day or several days, however, the Fermi- 
LAT data have been averaged over a three months pe- 
riod. Therefore, the multi-frequency data are quasi- 
simultaneous but not really simultaneous. Very recently, 
Giommi et al. (2011) built the SEDs for 105 blazars with 
the simultaneous Planck, Swift, Fermi, and ground-based 
data. The data of the 105 blazars are obtained between 
2009 Dec. and 2010 Oct., which were all collected in one 
day or several days, and the multi-frequency data are re- 
ally simultaneous. 

Since the Fermi data are publicly available, many au- 
thors have studied the physical properties of blazars with 
multi-frequency data from radio-X-ray to 7-ray bands. 
Ghisellini et al. (2010) constructed the SEDs for 85 blazars 
in LBAS, and found that there is a positive correlation 
between the jet power and the luminosity of the accre- 
tion disk in FSRQs. Later, Ghisellini et al. (2011) also 
analyzed the physical properties of blazars with high red- 
shifts in the 11 months Fermi AGN sample (the First LAT 
AGN Catalog: ILAC) by using the multi-frequency data. 
However, most of the multi-frequency data used by these 
authors are not simultaneous or quasi-simultaneous. 

In the study of 7-ray emissions of blazars, one of im- 
portant physical properties is the locations of the 7-ray 
emission regions. Many authors studied the locations of 
emission regions in different ways. For example, observed 
time lag of the 7-ray emission relative to broad emission 
lines for 3C 273 is used to constrain the 7-ray emission 
region Liu et al. (2011), the 7-ray absorptions by internal 
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photons of AGNs can be used to constrain the locations of 
the 7-ray emission regions (Sitarek & Bednarek 2008; Bai 
et al. 2009; Poutanen et al. 2010), and the 7-ray flux vari- 
ability of a blazar can infer its location of 7-ray emis- 
sion (Tavecchio et al. 2010). In this paper, we study the 
locations of the emission regions of 21 blazars whose si- 
multaneous and quasi-simultaneous multi-frequency SEDs 
have been given by Giommi et al. (2011) and Abdo et al. 
(2010a) in the frame of the multi-component model de- 
scribed in Dermer et al. (2009). 

In this study, wc use cosmological parameters 
{Ho,n^,flA) = (70 km s-i Mpc-\ 0.3, 0.7). 

2. Photon Emission and Absorption of FSRQs 

2.1. Photon Emission 

Broad band emission from blazars can be usually de- 
scribed in the framework of leptonic models. Since there 
is difficulty to explain SEDs of FRSQs in a simple one-zone 
SSC model (Abdo et al. 2010a), we study the photon emis- 
sion properties of FSRQs in the frame of one-zone, homo- 
geneous synchrotron and inverse Compton model given by 
Dermer et al. (2009). There are three soft photon seeds in 
this model: synchrotron photons, accretion disk photons, 
and broad line region (BLR) photons. 

Nonthermal photon emission is assumed to be produced 
by both the synchrotron radiation and IC scattering of 
relativistic electrons in a spherical blob of the jet which 
is moving relativistically at a small angle to our line of 
sight, and the observed radiation is strongly boosted by 
a relativistic Doppler factor Sd. The relativistic electron 
distribution which we used here is the same as that given 
by Dermer et al. (2009), i.e. 

7Vi(7') = Ki?(7';7;i„,7:nax){7'-"^ cxp(-777;) 
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where A'^^ is the normalization factor, H{x;xi,X2) is the 
Heaviside function: H{x;xi,X2) = 1 for a:i < x < a;2 and 
H{x;xi,X2) = everywhere else; as well as H{x) = for 
X < and H{x) = 1 for a; > 0. In the co-moving frame, this 
distribution is a double power law with two energy cutoffs: 
7^-i;jj and 7^nax- The spectrum is smoothly connected with 
indices pi and p2 below and above the electrons' break 
energy ^l^. Note that here and throughout the paper, un- 
primed quantities refer to the distant observer's frame on 
Earth and primed ones refer to the co-moving frame. 

After giving the electron distribution, the synchrotron 
flux can be given by Finke et al. (2008) 
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where e is the electron charge, B is the magnetic 
field strength, h is the Planck constant, and c^l is 
the distance to the source with a redshift z. Here 
TTieC^e' ~hv{l + z)l5-a is synchrotron photons energy 
in the co-moving frame, where mc is the rest mass 
of electron and c is the speed of light. In equation 



(2), R{x) = {x/2) j;;d0sm0 p^.^^dtK^/^dt, where x = 

ATTe'mlc^ /SeBh-f ^, 9 is the angle between magnetic field 
and velocity of high energy electrons, and Kc,^^[t) is the 
modified Bessel function of the second kind of order 5/3. 
Here we use an approximation for R(x) given by Finke 
et al. (2008). The synchrotron photons with vFj/^'^ have 
following spectral energy density 
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where 



is blob's radius and 
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variability timescale. 

For the IC scattering, there are two kinds of soft photon 
fields: non-thermal and thermal. The non-thermal field is 
the synchrotron photons field and high energy peak is cre- 
ated in the SSC process. For isotropic and homogeneous 
photon and electron distributions, the SSC flux, I'F^, is 
given by Finke et al. (2008) 
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where ctt is the Thomson cross section, mcC^e(, = hi/{l + 
2)/(5d is the energy of IC scattered photons in the co- 
moving frame, Fc{q' , T'J = 2q'lnq' + {1 + 2q'){l - q') + 
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The thermal soft photon fleld includes soft photon di- 
rect from accretion disk and broad line region (BLR) 
photon components and the IC process is called exter- 
nal Compton (EC) process. Here, we do not take the 
BLR line radiation into account but only consider the 
BLR Thomson scattered photon fleld. Assuming a stan- 
dard accretion disk with a central black hole mass M = 
A/glO^-Ajfsun and an accretion disk luminosity Lj, its mean 
photon energy in units of ra^f? at radius R is ed{R) = 



1.5 X 10-''( 
where 
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{F^^^e k Dermer 2010a), 



is the accretion efficiency, r, 



with LEdd = 1-26 X lO-^^Mgerg s"\ 77 
^ = 1.5 X lO^Mgcm, 



V3(i?) = l-/3i(-if )°■^ A = 1-0, and Rf^ = 6GM/c^ for the 
Schwarzschild metric. Therefore, the differential energy 
density per one solid angle of photon direct from accre- 
tion disk is 



tidisk(e,M*;''b) = K 



aid 



B7r(moC^) 



167r2?/(i?3/rg) 
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e«/ed(-R) _ 1 



where R = rbVu* ^ — 1, the normalization factor k = 
-T^^ , and Td = 1.164 x IQ'^taiRWcC^ is the tempera- 

ture of accretion disk at radius i?, and a is Boltzmann 
energy density constant (Dermer et al. 2009; Inouc & 
Takahara 1996). 
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For the BLR photon field, its specific energy density 
depends on the BLR structure and the radiation feature 
of accretion disk. Here, two assumptions are made: (1) 
the BLR is a spherically symmetric shell with inner ra- 
dius Ri and outer radius i?o- In general, the radius of 
BLR can be roughly estimated by using the relationship 
between some emission lines luminosity and the radius 
of BLR (e.g., Kaspi et al. 2005; Bentz et al. 2009), 
and the method of straight-forward reverberation map- 
ping (e.g., Paltani & Turler 2005), but as discussed in 
Ghiscllini & Tavecchio (2008), these methods also lead 
to large uncertainties on the radius of BLR. Fortunately, 
it has been found that the BLR radius i?BLR can scale 
roughly as L[J-^ (e.g., Wandel et al. 1999) or L^ '^ (e.g., 
Kaspi et al. 2000). Here, we assume the simplest hy- 
pothesis that the BLR radius scales with the square root 
of the disk luminosity (Ghisellini & Tavecchio 2008), i.e., 
-Rblr = 10^^[Ld/10''^ erg s~^]^/^ cm, consider the i?BLR as 
the outer radius Ro of the shell BLR, and assume ^ = C 
with ^ fa 10 representing a moderate width; and (2) the 
gas density of the BLR has the power-law distribution 
ne(r) = no(-^)^ for Ri < r < Ro, it has been suggested 
that the particle density of the BLR in quasars scales as 
J.-1.0 J.-1.5 (e,g,^ Kaspi & Netzer 1999), here we use 
C = — 1.0. Since the gas of BLR will also Thomson scatter 
the central accretion disk radiation to form the diffuse ra- 
diation from the BLR, the density no(r) can be estimated 

by using the radial Thomson depth tt = trx drric (r) 
(Dermer et al. 2009), where ctt is the Thomson cross sec- 
tion and r is the distance from the central black hole, here 
we use Tt = 0.01 which is the typical value (e.g., Donea 
& Protheroe 2003; Finke & Dermer 2010a). For the ra- 
diation field of standard accretion disk described by the 
multi-color blackbody spectrum, its luminosity in the sta- 
tionary frame is 



Lbb(e) = 87r^ 



2c(tooC^) 



X3 



£/ed(fl) _ 1 



RdR, 
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where Ac = h/rrtcC — 2.426 x 10" cm is the electron 
Compton wavelength and R^ut — SOOi?^^, is the outer ra- 
dius of the accretion disk. Therefore, the differential en- 
ergy density per one solid angle of BLR-scattered photon 
field is 



WBLR(e,/i*;?'b) = K 



-F{fi^,rh), 



(7) 
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where ro is classic electron radius, rb is the distance from 
the emission blob to the central black hole, and i^(/i, ,rb) is 
the function given by Dermer et al. (2009) (their Eq.(97)). 
From equation (7), the specific energy density of BLR- 
scattered photon field depends on the angle, 9^, between 
the directions of the photons and the electrons, where /Xh. = 
cos 6*,. 

The i^F^ spectrum of the EC process is given by 



^^EC(BLR,disk) ^ 2n 



(8) 



' , W(BLR,disk)(e>M*;?'b) 

de — X 



d'y 5 Ckorncl : 



+ (^)2. y = l-es/-f. 



' 7low ^ 

where Eg ~ hv /jtlcC^ and 7 = (5d7'- The inverse Compton 

kernel is Ckorncl = y + y-l - -L ^ J^^^2 

and e = 7e(l — ^*) (Dermer et al. 2009). 

For a FSRQ, the thermal emission may be significant, 
so we also take the thermal contribution ^/i^J^'"' from the 
standard accretion disk into account. Therefore, the local 
spectrum for total emission from a FSRQ is 

uFl°' = Z^Fjh- + [^F,«yn + j^FfC + ^ j;.EC(BLR,disk)] (g) 

2.2. Photon Absorption 

In above descriptions, we did not considered the ab- 
sorption of high energy photons. In fact, there are two 
possible absorption processes: internal and external. For 
the internal absorption, since the photons from BLR and 
accretion disk not only provide target radiation fields for 
IC process to produce 7-ray emission, but also attenuate 
7-rays through the pair-production process, the internal 
absorption should be taken into account. The internal ab- 
sorption optical depth in the radiation fields of the disk 
and the BLR. is given by Dermer et al. (2009) 



intl 



27r / dr 

'rb 



(1 — /i*)d^» 



(10) 



dea 



77 (^) ^ 



2/(E^(l+2;)(l-M.)) 

Mdisk(£, M* ; t) + mblr(£, ; r) 

where (T^^(s) is the 77 pair-production cross section (e.g., 
Gould & Schreder 1967; Dermer et al. 2009). On the 
other hand, there is also the absorption in the blob itself, 
i.e., the absorption due to the interaction with the inter- 
nal synchrotron radiation field (e.g. Sitarek & Bcdnarck 
2007). In this case, the internal absorption optical depth 
is given by Finke et al. (2008) 

[i + zfcndl 



-E^vFl^'^ 



(11) 



25 



2 777, c i V , m 1 11 ]]) 

1.236 X 10^" Hz and vF^^'^ is given by 



where v = 
Eq. 2. 

For the external absorption, because high energy pho- 
tons emitted from any FSRQ will interact with extragalac- 
tic background light (EBL) in the propagation in the ex- 
tragalactic space and then are absorbed through the pair- 
production process (e.g., Stecker et al. 1992; Finke et al. 
2010b), such an absorption should be taken into account. 
Here we adopt the EBL model of Finke et al. (2010b) to 
get the optical depth of a 7-ray photon with energy at 
redshift z, r™^ (£;.,, 2). 

Therefore the photons flux observed at the Earth is 

uF, = {vFl°')GM-r''''''{E^,z)~T"''\E^,r^)] (12) 



4 



Yan et al. 



Table 1. List of model parameters used to reproduce the SEDs of 21 FSRQs. The descriptions of columns in Tabic 1 are as follows: 
column[l]: name; column[2]: redshift; column[3]: magnetic field in Gauss; column[4-6]: minimum, break and maximum random 
Lorcntz factors of the rclativistic electrons; column[7 and 8]: slopes of relativistic electrons distribution below and above the 7^; 
column[9]: Doppler factor of the blob; column[10]: blob size in unit of cm; column[ll]: black hole mass in 10** times solar mass; 
column[12]: the accretion luminosity; column[13]: the distance from blob to black hole in cm. *: bad fits at optical or X-ray band, 
see text for the detailed comments. 



Name 


z 


B 


^min 


7; 


Tmax 


Pi 


P2 


Su 


K 


Ms 


^Edd 


rh 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[i 


i] 


[9] 


[10] 


[11] 


[12] 


[13] 


0133 + 47 


0.86 


0.15 


400 


3.2E3 


5.0E4 


2.0 


4, 


.0 


25 


8.13E16 


16.1 


0.03 


4.07E17 


0208-512* 


1.0 


0.18 


720 


3.3E3 


5.0E4 


2.0 


3 


.8 


24 


1.0E16 


8.29 


0.06 


3.25E17 


0227-369* 


2.12 


0.2 


980 


3.1E3 


5.0E4 


2.0 


4, 


.5 


25 


4.79E15 


10 


0.02 


2.03E17 


0347-211 


2.94 


0.22 


1.86E3 


3.5E3 


5.0E4 


2.0 


4, 


.8 


23 


2.1E16 


10 


0.06 


3.99E17 


0420 - 014* 


0.92 


0.18 


400 


3.0E3 


1.5E4 


2.0 


5, 


.0 


26 


1.23E16 


32.4 


0.01 


3.01E17 


0454 - 234* 


1.0 


0.2 


810 


3.1E3 


4.0E4 


2.0 


5, 


.0 


26 


4.37E15 


25.1 


0.01 


2.22E17 


0528 + 134 


2.07 


0.24 


1.3E3 


2.0E3 


2.2E4 


2.0 


5, 


.0 


23 


4.27E15 


45 


0.02 


8.99E17 


1454 - 345 


1.42 


0.24 


1.1E3 


2.oE3 


2.2E4 


2.0 


5, 


.0 


23 


1.58E16 


20 


0.02 


2.72E17 


4C 01.28 


0.89 


0.24 


380 


2.0E3 


9.0E3 


2.0 


4, 


.0 


23 


2.0E16 


17.8 


0.02 


2.84E17 


A 00 AT 

4C 28.07 


1.21 


0.24 


920 


2.0E3 


A ATP 

9.0E3 


A 

2.0 


4, 


.0 


23 


3.47E15 


1 A 

10 


A AO 

0.03 


2.29E17 


0917 + 449* 


2.19 


0.25 


760 


2.0E3 


9.0E3 


2.13 


4, 


.0 


32 


8.71E15 


75.8 


0.05 


9.81E17 


4C 29.45 


0.73 


0.22 


760 


2.0E3 


3.0E4 


2.13 


4, 


.0 


22 


4.57E15 


12.9 


0.02 


2.42E17 


3C 273 


0.16 


0.22 


210 


800 


1.25E3 


2.0 


5, 


.0 


16 


1.7E16 


10 


0.4 


7.6E17 


1510-089* 


0.36 


0.22 


355 


850 


1.7E3 


2.05 


5, 


.0 


16 


3.16E16 


13.5 


0.02 


1.7E17 


1308 + 32 


1.0 


0.18 


550 


3.3E3 


1.2E5 


2.0 


3 


.8 


24 


1.38E16 


3.98 


0.03 


1.61E17 


1502 + 106 


1.84 


0.18 


1.15E3 


3.3E3 


1.5E5 


2.0 


3 


.8 


26 


1.32E16 


31.6 


0.02 


3.13E17 


1520 + 319* 


1.49 


0.22 


355 


850 


3.7E3 


2.0 


5, 


.0 


16 


3.47E16 


50 


0.01 


1.38E17 


4C 66.20 


0.66 


0.18 


430 


3.3E3 


7.0E4 


2.0 


3 


.8 


26 


1.17E16 


41.5 


0.01 


3.09E17 


OX 169 


0.21 


0.18 


370 


3.3E3 


6.3E4 


2.0 


3 


.8 


26 


1.48E16 


9.54 


0.01 


2.2E17 


2325 + 093 


1.84 


0.15 


1.4E3 


2.0E3 


2.2E4 


2.0 


5, 


.0 


23 


5.37E15 


20 


0.03 


3.3E17 


2345- 1555* 


0.62 


0.15 


1.23E3 


2.0E3 


2.2E4 


2.0 


4, 


.0 


23 


2.14E15 


10 


0.01 


1.52E17 



^ / l-exp(-rint2) '\ 

where i^F*"* is given by Eq. (9). 
3. Applications 

Wc now apply the model describes in Section 2 to re- 
produce the quasi-simultaneous multi-frequency SEDs of 
21 FSRQs. The observed data of these 21 FSRQs are 
taken from Abdo et al. (2010a). In fact, Giommi et al. 
(2011) construct simultaneous multi-frequency SED of 105 
blazars, however, most of these blazars were observed for 
short time, and the numbers of collected high energy pho- 
tons are small. Hence, most of these blazars only have 
upper limit of the GeV flux, which are not helpful for our 
goal. Six FSRQs in the sample of Giommi et al. (2011) 
having three or four GeV data points are 0420-01, 0454- 
234, 4C 29.45, 3C 273, 1308+326 and 1502+106, which 
are also presented in the sample of Abdo et al. (2010a). 

The modelling results are shown in Figs. 1- 3. It can 
be found that our models can reproduce well the quasi- 
simultaneous emission from optical-X-ray to GeV band 
of most FSRQs. Our model can not reproduce very well 
the optical or X-ray spectrum of several sources (marked 
with * in Table 1). We believe that the bad reproduc- 
tion at optical band of PKS 1510-089 is caused by the 
fact that the observed optical data are contaminated by 



the emission from its host galaxy. For the source B2 
1520+31, it seems that the standard accretion disk model 
we used here does not adapt to describe its UV emission. 
The observed X-ray spectrum of these sources 0208-512, 
0227-369, 0420-014, 0454-234, 0917+449, and 2345-1555 
become flat at several keV band, which can not be re- 
produced by our model. In fact, the relativistic electron 
distribution with a large 7^iin can possibly explain such a 
spectrum (Tavecchio et al. 2009). The model parameters 
are listed in Table 1. The same values of?/ =1/12, C = — 1, 
tt = 0.01 and = 10 are used for all sources in our sam- 
ple, which are not listed in Table 1. The black hole masses 
of 0420-014, 1510-089 and 4C 66.20 (1849+64) are taken 
from Gu, Cao, & Jiang (2001). The black hole masses 
of 4C 01.28 (1055+018), 0917+449, 4C 29.45, 1308+32, 
1502-106 and OX 169 are taken from Chen et al. (2009). 
The black hole masses of the other sources are derived 
through modelling SEDs. It should be noted that the ther- 
mal emission of PKS 1510-089, 3C 273 and B2 1520+31 
are so significant that the optical-UV emission are purely 
from the thermal emission. Meanwhile, they have the 
smallest value of 7^, and J^^^- This indicates that the 
environment around the blob have a significant effect on 
the distribution of electrons in blob and the cooling of the 
electrons are dominated by EC process for these sources. 

From Table 1, we can find that the black hole masses 
are in the range (l-6)x lO^Afgun, the accretion luminos- 
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Fig. 1. Comparisons of modelling SEDs with the observed data for the Fermi FSRQs. For each FSRQ, thin solid line represents 
the thermal emission from accretion disk, dotted, dashed, and dash-dotted lines represent SSC, BLR-Compton, and disk-Compton 
components, respectively, and bold solid line is the sum of these components. The observed data are taken from Abdo et al. {2010a). 



ity of most of our sources are in the range 0.01 - 0.05 
(only one source 3C 273 requires the accretion luminosity 
larger than 0.1), and the disk luminosity cluster around 
6 X 10^^ erg s~^. For the high energy electron's distri- 
bution, we find that the values of the 7^;^ are in the 
range 400-1500, which are quite larger than the results 
derived by Ghisellini et al. (2010), indicating some clues 
to the acceleration and cooling processes of the electrons 
(Tavecchio et al. 2009; Katarzynski et al. 2006). From 
Table 1, it's found that 7max clusters around 3.0 x 10'', 
and besides three sources with ~ 800, 8 sources have 
7(j ~ 2000 and 10 sources have 7^ ~ 3000. For the magnetic 
field of radiation region, their values are distributed in a 
narrow range 0.15-0. 25 G, which are far smaller than the 
results found by Ghisellini et al. (2010). The small mag- 
netic field means inefficient synchrotron cooling, which are 
consistent with the results of large j'^^i^ and 7^, . The val- 
ues of Doppler factor 5u we derived range from 16 to 32, 
which is in the range reported by Savolainen et al. (2010). 

We now consider the relationships among size of the 
blob, I/d, and the location of the blob. It is expected 
that Tb^j relates to R'^ by = r^Oj, where 9^ is the open- 
ing angle of jet. However, we find that there is no corre- 
lation between R[^ and (see Fig. 4). The possible large 
distribution of 0j in our sample would disturb the expected 
correlation. On the other hand, this fact may support ei- 
ther the scenario of the 7-ray produced in a more compact 



region embedded in jet (Giannios, Uzdensky, & Begelman 
2009; Tavecchio et al. 2011) or the re-collimation of the 
flow (Bromberg & Levinson 2009). In our sample, we find 
that there is a strong correlation between and L<j with 
correlation coefficient r = 0.83 and a chance probability 
p = 1.94 X 10~^ (Fig. 5). This correlation is caused by the 
fact that rb cluster around 1.3i?o (top panel in Fig. 6). 

Our results confirm the well known cognitive that EC 
process is important for 7-ray emissions of the FSRQs. 
The contribution of EC component to 7-ray emission of 
the FSRQ is dominated by BLR component, which is de- 
termined by the location of the emission region (blob) 
rb and BLR structure. Therefore, 7-ray emission of the 
FSRQ can be used to constrain the location of emission 
blob. In Fig. 6, we present the distribution of the loca- 
tions of emission regions rb of 21 FSRQs. We find that 
all these sources have rb > 0.5i?o and the emission region 
of 19 out of 21 are outside of the BLR with the typical 
assumptions of BLR structure, i.e., their blobs are located 
in the distance range of 10^^ ^ 10^* cm. 

Since the EC component relates to the differential en- 
ergy density per one solid angle ublr of BLR photon field 
and wblr depends on the BLR structure. We now dis- 
cuss the uncertainty due to the parameters {(, tt, and 
Ro/Ri) describing the BLR structure. For these three pa- 
rameters, there is little observed information. Except for 
the typical values of BLR structure we used, there are al- 
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PKS 1454-354 



Fig. 2. Same as Fig. 1, but for different FSRQs. 





Fig. 3. Same as Fig. 1, but for different FSRQs. 
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15.0 15.5 



Fig. 4. The relationship between the location of the blob and 
the size of the blob. 



1^ 17.6 




45.0 45.5 46-0 46.5 47.0 47.5 
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Fig. 5. The relationship between the location of the blob and 
the accretion disk luminosity. 

tcrnativc values for these parameters. In a wind model 
for BLR suggested by Murray & Chiang (1995) and Elvis 
(2000), the value of ( should be -2. The alternative ex- 
treme values of tt is 0.1 (Dermer et al. 2009). The width 
of the BLR is poorly known, we will consider the cases of 
Ro/Ri = (2,10,60), which are corresponding to a very thin 
BLR shell, a moderate thin BLR shell and a thick BLR 
shell, respectively. From the observed 7-ray flux, the en- 
ergy density of BLR photon field required to reproduce 
the SEDs can be derived. 



.rh)de 



(13) 



We find that the values of C/BLR(?'b) are restricted in 
the range of (3.2 x lO"'^ - 7.9 x 10~'^)erg cm'"^ in our 
sample. Hence, no matter what assumption about the 
BLR structure we make, the allowed location of emission 
region would make the value of J7BLR(?'b) in the above 
range. In Fig. 7, we show the relationship between en- 




16.6 16.8 17.0 17.2 17.4 17.6 17.8 18.0 18.2 18.4 18.6 
Log r (cm) 



Fig. 6. Distributions of the locations of emission regions. 
Top panel: is in unit of outer radius of BLR Ro- Bottom 
panel: is in unit of cm. 

ergy densities of BLR-scattered photon field with differ- 
ent parameters for the BLR structure and the locations of 
emission regions with typical values in our sample, where 
we assume that A/rh = 2.0 x 10^ A^sun and i = 0.02 (i.e. 
id ~ 5 X 10"*^ erg s^^). We find that the smallest allowed 
value of Tb « 7.9 x 10^^ cm is derived when tt = 0.01, 
C = -2 and i?o/i?i = 60. With tt = 0.1, C = (-1, -2) and 
Ro/Ri = 2, the largest allowed value of rb ~ 1.3 x 10^^ 
cm is obtained. Since there are three sources with high 



accretion disk luminosity (id > 10*^ erg s~^) in our sam 
pie, and the accretion disk luminosity has a affect on con- 
straining the location of blob, we also discuss the case of 
high accretion disk luminosity Ld ~ 6.3 x 10"*^ erg s^^ with 
Mbh = 5.0 X 10^ Afsun and ^ = 0.1. The resuhs are shown 
in Fig. 8. It is found that the allowed value of rb is in the 
range 2.6 x 10^^ ~ 4.2 x 10^® cm. 

Finally, wc consider the roles of both internal and ex- 
ternal absorptions. For the internal absorptions in our 
sample, we find that they are negligible in reproducing 
the SEDs. As examples, in Fig. 9, we show the inter- 
nal optical depths in the radiation fields of the disk and 
the BLR for three sources (B2 1520-1-31, PKS 0528-1-134, 
and OX 169) and those in the synchrotron radiation fields 
for PKS 0528-^134, PKS 0347-211, and PKS 0227-369. 
From Fig. 9, B2 1520-f31 has a largest value of r'"" 
due to its high energy density of external photon field, 
^blr(''™'") ~ '^■^ ^ 10^"^ erg cm""^, which becomes signif- 
icant for photons with ly > 10^^ Hz. As to PKS 0528-1-134 
and OX 169 in Fig. 9, the former has the intermediate 
optical depth and the latter has the lowest optical depth 
in our sample. We expect that the 7-ray spectrum be- 
tween 0.1 - 10 GeV could not be reproduced well if the 
energy density of BLR photon field exceeds U^£^{r^^"^). 
For the absorption in blob itself, we show largest values 
of the optical depths, r'"*^, in our sample, it can be seen 
that these absorptions are negligible in our sample. On 
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Fig. 7. Change of energy density of BLR-seattered photon 
field with the location ri, of emission region for different as- 
sumptions of the BLR structure. The black points are the 
value obtained from the SED fits of modelled sources. The 
horizonal lines are the lower and upper limits of !7blr de- 
rived from the 7-ray emission of our sample, respectively. 
Dashed, dotted, and dash-dotted lines correspond to the cases 
of Ro/R{ = 2 , 10, and 60, respectively. For accretion lumi- 
nosity and black hole mass, the typical values found in our 
sample arc used: A/bh = 2.0 X 10*' A/sun and t = 0.02. 

the other hand, the fact that there is no significant spec- 
tral cutoff between f - f 00 GeV for FSRQs (Costamante 
2011; Ackerman et al. 2011) also supports the minimum 
value of 7'b we derived. For the external absorption, it de- 
pends on the EBL models, so different models will make 
some differences, in particular, for higher redshifts. We 
found that 15 sources have the maximum gamma-ray en- 
ergy i?niax ^ 6.6 GeV, 5 sources have £^max ~ 20.8 GeV, and 
1 source has ii'max ~ 65.7 GeV in our sample. Therefore, 
we calculate the external optical depths for six sources 
(4C 66.20, 4C 01.28, PKS 0454-234, PKS 1454-354, PKS 
1502-f 106, and S4 0917-^44 ) with E^^^ > 20 GeV us- 
ing four kinds of the EBL models given by Kneiske et al. 
(2004), Stecker et al. (2006), Franceschini et al. (2008), 
and Finke et al. (2010), respectively, and show the results 
in Fig. 10. It can be seen that our modelled SEDs for three 
sources with z < 1 are not affected by the EBL models, but 
the SEDs for the sources (PKS 1454-354, PKS 1502+106, 
and S4 0917-1-44 ) with z > 1.4 will be different for dif- 
ferent EBL models, particularly, for the EBL model given 
by Stecker et al. (2006). However, the EBL model given 
by Stecker et al. (2006) has been ruled out by Fermi-LAT 
observations (Abdo et al. 2010c). For PKS 1502-^106, 
the optical depth given by Kneiske et al. (2004) could be 
slightly larger than that expected by Finke model, which 
would lead to slightly difference for modeling the SED 
of PKS 1502-F-106. Considering the specific SED of PKS 
1502-1-106, we expect that this difference has no effect on 
our results. 



Fig. 8. Change of energy density of BLR-scattered photon 
field with the location of emission region for different as- 
sumptions of the BLR structure. The symbols are the same 
as that in Fig. 7. For accretion luminosity and black hole 
mass, Mbh = 5.0 X 10^ Msun and £ = 0.1 are used. 




v(Hz) 



Fig. 9. The optical depths of internal absorption in modeling 
the SEDs of FSRQs. Solid line: the absorption in the accre- 
tion disk and BLR radiation fields t'°*^ (from left to right: 
B2 1502-1-106, PKS 0528-1-134, and OX 169, which has the 
largest, intermediate and the lowest optical depth of r'"*^, 
respectively. B2 1502-1-106, the source with the strongest ab- 
sorption in the accretion disk and BLR radiation fields, is 
marked with a thick curve.); Dashed line: the absorption in 
blob ri'i'^ (fj.ojjj Ycft to right: PKS 0528-1-134, PKS 0347-211, 
and PKS 0227-369, which have the largest optical depths of 
^int2 y discussion in text. 



4. Conclusions 

In this paper, we have modeled the quasi-simultaneous 
multi-band spectra of 21 Fermi FSRQs in the frame of a 
multi-component one-zone leptonic emission model, and 
studied the locations of emission regions of the FSRQs 
in particular. We have found that the emission regions 
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Fig. 10. Change of the photon energy at T^^^{E-y,z) = 1 
with redshift z for different EBL models. Solid, dashed, 
dotted, and dot-dashed lines represent the results by using 
the models given by Finke et al. (2010), Franceschini et al. 
(2008), Knoiske et al. (2004), and Stocker et al. (2006). Filled 
diamonds represent the maximum 7-ray photon energy for six 
sources (4C 66.20, 4C 01.28, PKS 0454-234, PKS 1454-354, 
PKS 1502-1-106, and S4 0917-1-44 ). 

lie in the region of 7.9 x lO^^ - 1.3 x 10^** cm (300 - 
4300 Schwarzschild radii) for FSRQs with low accretion 
disk luminosity, and the emission regions locate in the 
larger region of 2.6 x 10^^ - 4.2 x 10^^ cm (300 - 5600 
Schwarzschild radii) for ones with high accretion disk lu- 
minosity. Our results are consistent with the results of 
Tavecchio et al. (2010) derived by analyzing the variabil- 
ity timescale and that of Liu, Bai, & Ma (2008) and Bai 
et al. (2009) obtained by discussing internal absorption. 
Our results disfavor the far dissipation scenario suggested 
by e.g., Marscher et al. (2010), in which 7-rays are pro- 
duced at larger distance (>10 pc) from black hole. At 
such large distance, the external radiation field for IC pro- 
cess are dominated by the emission from the dusty torus 
with fiR = 3.0 X 10^3 Hz and J7dusty < 3 x 10"^ erg cm'^ 
(Ghisellini & Tavecchio 2008). As we discussed above, 
the observed 7-ray flux requires that the energy density 
of external radiation field for IC process should be larger 
than 3.2 x 10~* erg cm~'^ with the moderate Doppler fac- 
tor, which is inconsistent with the far dissipation scenario. 
Moreover, the observed variability on timescale of ^ 10 
min of PKS 1222-f 216 (Aleksic et al. 2011; Tavecchio et 
al. 2011) requires quite fast IC cooling, which also chal- 
lenges the far dissipation scenario considering lower energy 
infrared photons as seed photons. 

We did not consider the target photons from BLR emis- 
sion line in this work. Luminosity of the BLR emission 
line is comparable with or larger than that of the BLR 
Thomson scatted photons radiation between 10^^- 10^^ Hz 
in the BLR zone (Tavecchio & Ghisellini 2008; Sbarrato 
et al. 2011), which would make the locations of the emis- 
sion regions we derived closer to the central black hole 
than the real ones. Based on the above discussion, we can 
roughly estimate that how much the emission line radia- 



tion affect our results. From our results, we can see that 
when the energy density of BLR-scattered photons field 
increases by '--^lO times (corresponding to varying from 
5 X 10^5 to 6.3 X 10^6 erg s-\ Figs. 7 and 8), the range 
of the emission region locations increases by ^^3 times. 
Therefore, if the energy density of the BLR emission line 
photon field is ~10 times that of the BLR Thomson scat- 
ted photons field, we expect the locations of the emission 
regions possibly will locate in the range 3.0 • (7.9 x 10^^ - 
1.3 X 10^*) cm and 3.0 • (2.6 x lO" - 4.2 x 10^^) cm for low 
and high accretion luminosity FSRQs, respectively. The 
detail contribution of BLR emission line will be studied 
in the future work. 
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